Non-minimal curvature-matter couplings in modified gravity 
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Recently, in the context of f{R) modified theories of gravity, it was shown that a curvature- 
matter coupling induces a non-vanishing covariant derivative of the energy-momentum, implying 
non-geodesic motion and, under appropriate conditions, leading to the appearance of an extra force. 
We study the implications of this proposal and discuss some directions for future research. 
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I. INTRODUCTION 



Current experimental evidence indicates that gravita- 
tional physics is in agreement with Einstein's theory of 
General Relativity (GR) to considerable accuracy (for 
thorough discussions see [1]); however, quite fundamen- 
tal questions suggest that it is unlikely that GR stands 
as the ultimate description of gravity. Actually, difhcul- 
ties arise from various corners, most particularly in con- 
nection to the strong gravitational field regime and the 
existence of spacetime singularities. Quantization is a 
possible way to circumvent these problems, nevertheless, 
despite the success of gauge field theories in describing 
the electromagnetic, weak, and strong interactions, the 
description of gravity at the quantum level is still miss- 
ing, despite outstanding progress achieved, for instance, 
in the context of superstring/M-theory. 

Furthermore, in fundamental theories that attempt to 
include gravity, new long-range forces often arise in ad- 
dition to the Newtonian inverse-square law. Even if one 
assumes the validity of the Equivalence Principle, Ein- 
stein's theory does not provide the most general way to 
establish the spacetime metric. There are also impor- 
tant reasons to consider additional fields, especially scalar 
fields. Although the latter appear in unification theories, 
their inclusion predicts a non-Einsteinian behaviour of 
gravitating systems. These deviations from GR include 
violations of the Equivalence Principle, modification of 
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large-scale gravitational phenomena, and variation of the 
fundamental couplings. 

On large scales, recent cosmological observations lead 
one to conclude that our understanding of the origin and 
evolution of the Universe based on GR requires that 
most of the energy content of the Universe is in the 
form of currently unknown dark matter and dark energy 
components that may permeate much, if not all space- 
time. Indeed, recent Cosmic Microwave Background 
Radiation (CMBR) data indicate that our Universe is 
well described, within the framework of GR, by a nearly 
flat Robertson- Walker metric. Moreover, combination of 
CMBR, supernovae, baryon acoustic oscillation and large 
scale structure data are consistent with each other only if, 
in the cosmic budget of energy, dark energy corresponds 
to about 73% of the critical density, while dark matter to 
about 23% and baryonic matter to only about 4%. Sev- 
eral models have been suggested to address issues related 
to these new dark states. For dark energy, one usually 
considers the so-called "quintessence" models, which in- 
volves the slow-roll down of a scalar field along a smooth 
potential, thus inducing the observed accelerated expan- 
sion (see [2] for a review). For dark matter, several 
weak-interacting particles (WIMPs) have been suggested, 
many arising from extensions to the Standard Model {e.g. 
axions, neutralinos). A scalar field can also account for 
an unified model of dark energy and dark matter [3] . Al- 
ternatively, one can implement this unification through 
an exotic equation of state, such as the generalized Chap- 
lygin gas [4]. 

However, recently a different approach has attracted 
some attention, namely the one where one considers a 
generalization of the action functional. The most straight 
forward approach consists in replacing the linear scalar 
curvature term in the Einstein-Hilbert action by a func- 
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tion of the scalar curvature, f(R). In this context, a 
renaissance of f{R) modified theories of gravity has re- 
cently been verified in an attempt to explain the late- 
time accelerated expansion of the Universe (see for in- 
stance Refs. [5, 6] for recent reviews). One could alter- 
natively, resort to other scalar invariants of the theory 
and necessarily analyze the observational signatures and 
the parameterized post-Newtonian (PPN) metric coeffi- 
cients arising from these extensions of GR. In the con- 
text of dark matter, the possibility that the galactic dy- 
namics of massive test particles may be understood with- 
out the need for dark matter was also considered in the 
framework of f{R) gravity models [7]. Despite the ex- 
tensive literature on these f{R) models, an interesting 
possibility has passed unnoticed till quite recently. It 
includes not only a non-minimal scalar curvature term 
in the Einstein-Hilbert Lagrangian density, but also a 
non-minimal coupling between the scalar curvature and 
the matter Lagrangian density [8] (see also Ref. [9] for 
related discussions) . It is interesting to note that nonlin- 
ear couplings of matter with gravity were analyzed in the 
context of the accelerated expansion of the Universe [10], 
and in the study of the cosmological constant problem 
[11]. In this contribution we discuss various aspects of 
this proposal. 

This work is organized as follows: in the following Sec- 
tion, the main features of this novel model are presented. 
In Section III, the issue of the degeneracy of Lagrangian 
densities, actually a feature well known in GR [12-14], is 
addressed in the context of the new non-minimally cou- 
pled model [15]. In Section IV, the scalar-tensor repre- 
sentation of the model is presented, with particular em- 
phasis on the new features and difficulties encountered in 
the new model. These issues arc quite relevant, as they 
allow one to properly obtain the PPN parameters (3 and 
7 and show that they arc consistent with the observations 
[16]. Section V, addresses the compatibility of the model 
with the astrophysical condition for stellar equilibrium 
[17]. In Section VI, a further generalization of the model 
is discussed and an upper bound on the extra acceleration 
introduced by the new non-minimal coupling is obtained 
[18]. Finally, in Section VI our conclusions are presented 
and objectives for further research are discussed. 

Throughout this work, the convention SttG = 1 and 
the metric signature (— , -|-, +, +) are used. 



II. LINEAR CURVATURE-MATTER 
COUPLINGS 

The action for curvature-matter couplings, in f{R) 
modified theories of gravity [8] , takes the following form 



-I 



-h{R) + [l + Xf2{R)]jCr. 



(1) 



where fi{R) (with i = 1,2) are arbitrary functions of 

the curvature scalar R and Cm is the Lagrangian density 
corresponding to matter and A is a constant. Since the 



matter Lagrangian is not modified in the total action, 
these may be called modified gravity models with a non- 
minimal coupling between matter and geometry. 

Varying the action with respect to the metric g^v yields 
the field equations, given by 



FiR^,,. - - V^V.Fi + g^^DF^ = (1 + A/2)T^. 



-2\F2CmRnv + 2A(V^Vi/ — gij.i'^)CmF2 ■ 



(2) 



where one denotes Fi{R) = f[{R)., and the prime denotes 
differentiation with respect to the scalar curvature. The 
matter energy-momentum tensor is defined as 



T 



2 6(^Cm) 
/=5 Sigi^'^) 



(3) 



Now, taking into account the generalized Bianchi identi- 
ties, one deduces the following generalized covariant con- 
servation equation 



l + A/2 



[9i,.jCm - T^.] VR. 



(4) 



It is clear that the non-minimal coupling between curva- 
ture and matter yields a non-trivial exchange of energy 
and momentum between the geometry and matter fields 
[16]. 

Considering, for instance, the energy-momentum ten- 
sor for a perfect fluid. 



Ti^u = {p + p) U^Uv + pg^v , 



(5) 



where p is the energy density and p is the pressure, re- 
spectively. The four- velocity, 17^, satisfies the conditions 
UfjJJ^ = —1 and U^Un-v = 0. Introducing the projec- 
tion operator h^i, = g^i, -\- U^U^, one can show that the 
motion is non-geodesic, and governed by the following 
equation of motion for a fluid element 



ds ' '"^ 

where the extra force, appears and is given by 



(6) 



p + p 



XFo 



l + A/i 



{Cm - p) V^i? + V^p 



h'"' . (7) 



One verifies that the first term vanishes for the specific 
choice of Cm — P, as noted in [19]. However, as pointed 
out in [15], this is not the unique choice for the La- 
grangian density of a perfect fluid, as will be outlined 
below. 



III. PERFECT FLUID LAGRANGIAN 
DESCRIPTION 

The novel coupling in action (1) has attracted some 

attention and, in a recent paper [19], this possibility has 
been applied to distinct matter contents. Regarding the 
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latter, it was argued that a "natural choice" for the mat- 
ter Lagrangian density for perfect fluids is = P, based 
on [12, 13], where p is the pressure. This specific choice 
implies the vanishing of the extra force. However, al- 
though £m = P does indeed reproduce the perfect fluid 
equation of state, it is not unique: other choices include, 
for instance, Cm = —p [13, 14], where p is the energy 
density, or = —na, where n is the particle num- 
ber density, and a is the physical free energy defined as 
a = p/n — Ts, with T being the fluid temperature and s 
the entropy per particle. 

In this section, following [13, 15], the Lagrangian for- 
mulation of a perfect fluid in the context of GR is re- 
viewed. The action is presented in terms of Lagrange 
multipliers along the Lagrange coordinates a"* in order 
to enforce specific constraints, and is given by 

Sm = j d'^x [-^/^ p{n, s) + J" (tp,^ + s^,^ + /^Aaj,)] • 

(8) 

Note that the action Sm = S{gni,, J^, (p, 9, s, a^, (3a) is a 
functional of the spacetime metric ^j^y, the entropy per 
particle s, the Lagrangian coordinates a^, and spacetime 
scalars denoted by ip, 6, and (3 a, where the index A takes 
the values 1, 2, 3 (see [13] for details). 

The vector density J'' is interpreted as the flux vec- 
tor of the particle number density, and defined as J** = 
^—gnlJ^. The particle number density is given by 
n = \J\/^/^, so that the energy density is a function 
p = p{\J\/\/^, s). The scalar field ip is interpreted as a 
potential for the chemical free energy /, and is a Lagrange 
multiplier for J'^, the particle number conservation. The 
scalar fields are interpreted as the Lagrange multipli- 
ers for a'^ = 0, restricting the fiuid 4— velocity to be 
directed along the flow lines of constant a^. 

The variation of the action with respect to J**, ip, 6, 
s, and (3a, provides the equations of motion, which 
are not written here (we refer the reader to Ref. [15] for 
details). Varying the action with respect to the metric, 
and using the definition given by Eq. (3), provides the 
stress-energy tensor for a perfect fluid 



dp 
dn 



with the pressure deflned as 



dp 

P = n- p. 

on 



p (5'"' + V'U-) , (9) 



(10) 



This definition of pressure is in agreement with the First 
Law of Thermodynamics, dp = pdn + nTds. The lat- 
ter shows that the equation of state can be specified by 
the energy density p(n,s), written as a function of the 
number density and entropy per particle. The quantity 
H = dp/dn = {p + p)/n is defined as the chemical po- 
tential, which is the energy gained by the system per 
particle injected into the fluid, maintaining a constant 
sample volume and entropy per particle s. 



Taking into account the equations of motions and the 
definitions ~ \f-9 nU^ and p. = {p + p) /n, the ac- 
tion Eq. (8) reduces to the on-shell Lagrangian density 



r 



: p, with the action given by [15] 



Sm — I d X 



(11) 



which is the form considered in Ref. [12]. It was a La- 
grangian density given by Cm = P that the authors of 
[19] use to obtain a vanishing extra-force due to the non- 
trivial coupling of matter to the scalar curvature R. For 
concreteness, replacing Cm = p in Eq. (7), one arrives at 
the general relativistic expression 



r 



P + P 



(12) 



However, an on-shell degeneracy of the Lagrangian 
densities arises from adding up surface integrals to the 
action. For instance, consider the following surface inte- 
grals added to the action Eq. (8), 

-/rf4ar(<^J''),^ , - jd^xiOsJ'')^^, 
-Jd^x{J'^(3Aa^),^ , 

so that the resulting action takes the form 



-gp{n, s) - ipJl"^ 



- e{sJ^),^-a^{PAJn,^^ 



(13) 



This action reproduces the equations of motion, and tak- 
ing into account the latter, the action reduces to [15] 



Sm — I d X 



(14) 



i.e., the on-shell matter Lagrangian density takes the 
following form Cm = —p- This choice is also consid- 
ered for isentropic fluids, where the entropy per parti- 
cle is constant s = const. [13, 14]. For the latter, the 
First Law of Thermodynamics indicates that isentropic 
fiuids are described by an equation of state of the form 
a{n, T) = p{n)/n-sT [13] (see Ref. [20] for a bulk-brane 
discussion of this choice). 

For this specific choice of Cm{2) = —p the extra force 
takes the following form: 



l + A/2 



1 



■P 



(15) 



An interesting feature of Eq. (15) is that the term related 
to the specific curvature-matter coupling is independent 
of the energy-matter distribution. 

The above discussion confirms that if one adopts a par- 
ticular on-shell Lagrangian density as a suitable func- 
tional for describing a perfect fluid, then this leads to 
the issue of distinguishing between different predictions 
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for the extra force. It is therefore clear that no straight- 
forward conclusion may be extracted regarding the ad- 
ditional force imposed by the non-minimal coupling of 
curvature to matter, given the different available choices 
for the Lagrangian density. One could even doubt the 
validity of a conclusion that allows for different physi- 
cal predictions arising from these apparently equivalent 
Lagrangian densities. 

Despite the fact that the above Lagrangian densities 
jC^(^i) are indeed obtainable from the original action, it 
turns out that they are not equivalent to the original La- 
grangian density Cm- Indeed, this equivalence demands 
that not only the equations of motion of the fields de- 
scribing the perfect fluid remain invariant, but also that 
the gravitational field equations do not change. Indeed, 
the guiding principle behind the proposal first put for- 
ward in Ref. [8] is to allow for a non-minimal coupling 
between curvature and matter. 

The modification of the perfect fluid action Eq. (8) 
should only affect the terms that show a minimal cou- 
pling between curvature and matter, i.e., those multi- 
plied by [15]. Thus, the current density term, which 
is not coupled to curvature, should not be altered. Writ- 
ing Cc = -p{n, s),Vu, = (p^^ + s9^^+(3Aaf^, for simplicity, 
the modified action reads 

^[l + Xf2{R)]Ce + J''V^ + Bll] , 

(16) 

and one can see that only the non-minimal coupled term 
Cc appears in the field equations, as variations with re- 
spect to g^^"^ of the remaining terms vanish: 

FiR^, - ^hgf,, - V^V.Fi + g^.DFi = (1 + A/2)T^, 
-2XF2CcR^. + 2A(V^V, - g^,D)CcF2 . (17) 

Thus, quite logically, one finds that different predic- 
tions for non-geodesic motion are due to different forms 
of the gravitational field equations. Therefore, the equiv- 
alence between different on-shell Lagrangian densities 
Cm(i) and the original quantity Cm is broken, so that one 
can no longer freely choose between the available forms. 
For the same reason, the additional extra force is unique, 
and obtained by replacing Cc = —p into Eq. (7), yielding 
expression (15). 

Indeed, in a recent paper [29], a generalization of 
the above approach is considered, by using a systematic 
method that is not tied up to a specific choice of matter 
Lagrangians. In particular, the propagation equations 
for pole-dipole particles for a gravity theory with a very 
general coupling between the curvature scalar and the 
matter fields is examined, and it is shown that, in gen- 
eral, the extra-force does not vanish. 



IV. SCALAR-TENSOR REPRESENTATION 

The connection between f(R) theories of gravity and 
scalar-tensor models with a "physical" metric coupled to 



the scalar field is well known. In this section, one pursues 
the equivalence between the model described by Eq. (1) 
and an adequate scalar-tensor theory. In close analogy 
with the equivalence of standard f{R) models [21], this 
equivalence allows for the calculation of the PPN param- 
eters (3 and 7 [22]. 

One may first approach this equivalence by introducing 
two auxiliary scalars ip and </> [19], and considering the 
following action 



Si 



I 



i/l(0) + [1 + A/2(0)] Cm + - <!>) 



(18) 

Now, varying the action with respect to V gives (f) = R 
and, consequently, action (1) is recovered. Varying the 
action with respect to </>, yields 



i) = -Fi+ XF2Cr, 



(19) 



Substituting this relationship back in (18), and assuming 
that at least one of the functions fi is nonlinear in R, one 
arrives at the following modified action 



-I 



[1 + A/2 (</.)]£„ 



-Fi{<t>) + \F2{<t>)C. 



-g rf'a(?o) 



where one still verifies the presence of the curvature- 
matter coupling. Note that this is not an ordinary scalar- 
tensor theory, due to the presence of the third and last 
terms. The former represents a scalar-matter coupling, 
and the latter a novel scalar-curvature-matter coupling. 
One may also use alternative field definitions to cast the 
action (18) into a Bran-Dicke theory with w = 0, i.e. no 
kinetic energy term for the scalar field, but with the ad- 
dition of a i?-matter coupling [19]. In conclusion, despite 
the fact that the introduction of the scalar fields helps in 
avoiding the presence of the nonlinear functions of R, the 
curvature-matter couplings are still present and, conse- 
quently, these actions cannot be cast into the form of a 
familiar scalar-tensor gravity [19]. 

However, one may instead pursue an equivalence with 
a theory with not just one, but two scalar fields [16]. 
This is physically well motivated, since the non-minimal 
coupling of matter and geometry embodied in Eq. (1) 
gives rise to an extra degree of freedom (notice that the 
case of a minimal coupling /2 = yields -(/i = Fi((/))/2, so 
that this degree of freedom is lost). Indeed, action of Eq. 
(18) may be rewritten as a Jordan-Brans-Dicke theory 
with a suitable potential, 



Si 



I 



VJ?-y(<^,^^) + [l + A/2((^)]£. 



-g d^x , 



with !/(</., V) =</>V'-/i('/')/2. 



(21) 
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Variation of this action yields the field equations 



V' ' 



(22) 



1 VU,i;) 1 „ , 



V' 



which, after the substitutions 4> = R and = Fi/2 + 
\F2Cm, collapses back to Eqs. (2). Likewise, the Bianchi 
identities yield the generalized covariant conservation 
equation 



VT^. = - 



A/2 



(0 - R) V^V + 



also equivalent to Eq. (4). 

Through a conformal transformation g^i, — > g*^^^ = 
tj^fjfiu (sec e.g. [23]), the scalar curvature can decouple 
from the scalar fields, so that the action is written in the 
so-called Einstein frame). A further redefinition of the 
scalar fields, 



1 V3, , 



^^=^, (24) 
allows the theory to be written canonically, that is. 



Si 



-I 



R* - 23*^ Vi,-^:^^, 



(25) 



V-g* d'^x, 



with = jCrn/i'^, the redefined potential 



2V3 1 



f X 



(26) 



and the metric in the field space {(fi^,ip'^), 



1 





(27) 



which, after a suitable addition of an anti-symmetric 
part, will be used to raise and lower Latin indexes. 

Variation of action Eq. (25) with respect to the metric 
dfiv yields the field equations 

K'^ - \9luR* = SttG (1 + A/2) T;, + (28) 



while variation with respect to gives the Euler- 

Lagrange equations for each field: 

a*(p' =B'+ AttG [a' (1 + A/2) T* - Xct'^F2£*] (29) 
where one defines Bi = dU /d(p^ and 

"^ = ~2^^^ ^ "^=-^ ' "^=°' (^^^ 

Eqs. (28), together with the Bianchi identities, result 
in the generalized conservation law 



(23) V*^T, 



AFo 



1 + A/2 



(31) 

Prom current bounds on the Equivalence Principle, it is 
reasonable to assume that the effect of the non-minimum 

coupling of curvature to matter is weak, A/2 <C 1. Sub- 
stituting this into (31) one gets, at zeroth-order in A, 



(32) 



so that one may disregard the /2(9?^) factor in the action 
(25) and consider only through the coupling present in T* 
(stemming from the definition of C^) and the derivative 
of tf^ (since cx \ogip and ij} = Fi + F2C). 

If both scalar fields are light, leading to long range 
interactions, one may calculate the PPN parameters (3 
and 7 [22], given by 



a a-' a 



(l + a2 



7-1 





(33) 

where aj^i = daj/dip'^ and = a^a'^ = a'^^a^aj; the 
subscript refers to the asymptotic value of the related 
quantities, which is connected to the cosmological values 
of the curvature and matter Lagrangian density. Prom 
the values found in Eq. (30), one concludes that /3 = 7 = 
1, as obtained in GR. However, it should be expected that 
small deviations of order 0(A) arise when one considers 
the full impact of Eq. (31). 

Furthermore, it should be empathized that the added 
degree of freedom embodied in the non-minimal /2 ^ 
coupling is paramount in obtaining values for the PPN 
parameters f3 and 7 within the current experimental 
bounds (or, conversely, allowing for future constraints of 
the magnitude of A and the form of /2); indeed, in the 
case where only the curvature term is non-trivial, fi^R 
and /2 = 0, one degree of freedom is lost and the param- 
eter a j^O defined in Eq. (30) is no longer a vector, but a 
scalar quantity: as a result, 7^ and one gets 7 = 1/2. 
In the discussed model, the vector aj has = 0, thus 
solving this pathology (see [16] and references therein for 
a thorough discussion) . 

Finally, notice that these results should be indepen- 
dent of the particular scheme chosen for the equivalence 
between the original model and a scalar-tensor theory; 
this may be clearly seen by opting for a more "natural" 
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choice for the two scalar fields (in the Jordan frame), 
such that (p — R and %!) = £. Although more physically 
motivated, this choice of fields is less pedagogical and 
mathematically more taxing [16]. 



V. IMPLICATIONS FOR STELLAR 
EQUILIBRIUM 

In this section, one studies the impact of the non- 
minimally coupled gravity model embodied in action Eq. 

(1) in what may be viewed as its natural proving ground: 
regions where curvature effects may be high enough, 
to evidence some deviation from GR, although moder- 
ate enough so these are still perturbative - a star [17] 
(see also [24] for other physical examples of the adopted 
methodology). As will be shown, the purpose of this 
exercise is to calculate deviations to the central temper- 
ature of the Sun (known with an accuracy of 6%), due 
to the perturbative effect of the non-minimal coupling of 
geometry to matter. 

Clearly, a full treatment of the equations of motion 

(2) is exceedingly demanding, unless a specific form for 
fi{R) and f2(R) is considered. Furthermore, since one 
is mainly interested in the ascertaining the effects of the 
non-minimal coupling within a high curvature and pres- 
sure medium, the modifications due to the pure curvature 
term /i should be overwhelmed by the effect of /2; under 
such circumstances, one may discard the former term, as 
thus take the trivial f\ = R case. A thorough discussion 
on the validity of this approximation with regard to rep- 
resentative, physically viable candidates for the function 
fi{R) is found in Ref. [17]. 

One now deals with the particular form of the cou- 
pling function /2. One considers the simplest form, which 
might arise from the first order expansion of a more gen- 
eral function in the weak field environ of the Sun, f2 = R 
(this implies that [A] — M~^). Also, one assumes that 
stellar matter is described by an ideal fluid characterized 
by a Lagrangian density Cm = P, [12, 13]. Adopting 
/i = /2 = R, the field equations become 

(1 + 2Xp) R^, - ^R {g^, + 2AV) = (34) 

2A(V^V. -5^.n)p + 

Notice that both \p and Xp are dimensionless quantities: 
the perturbative condition A/2 <C 1 translates to Ap <C 1 
and Xp <C 1. 

Taking the trace of the above equation yields 



R 



'Sp — p + (iAHp 
2[l + A(p-5p)] 



(35) 



inserting T = Tj^ = p — 3p. Substituting this into Eq. 
(34) and keeping only first order terms in A, one obtains 



2[H-A(p-3p)]i?^, = (36) 
(3p - p)gi.u + 2(1 - 2A?>) V + 2A(4V^V. - gi,.U)p , 



Since temporal variations are assumed to occur at the 
cosmological scale H^^ ^ , and are thus negligible at an as- 
trophysical time scale, one considers an ideal, spherically 
symmetric system, with a line element derived from the 
Birkhoff metric (in its anisotropic form) 



ds 



sm 



(e<^('-)dr2+df}2j , (37) 
'9 d(p'^). Following the usual 



with dVl = 

treatment, one defines the effective mass rUe through 
e~'^ = 1 — 2Gme/r which, replacing in Eq. (36) yields, 
to first order in A, 



m'^^A^r^p [l + 2Afp-^-|^ 

L V 2 2 p 

^ ( 5e-"VoVo + ^e-^VrVr + 2 



+ 



(38) 



VgVe \ 
r2 ) 



P 



which clearly shows the perturbation to the gravitational 
mass, defined by m'g = Anr'^p (in here, the prime denotes 
differentiation with respect to r). 
Taking the Newtonian limit 



> 2Gm 



;(r) , p{r) ^ p(r) , me(r) ^ 47rp(r)r'^ 



(39) 



and going through a few algebraic steps (depicted in [17]), 
one eventually obtains the non-relativistic hydrostatic 
equilibrium equation 



P 



•' + 



GrUe 



2A 



^P 



AwGpp 



p + PP 



(40) 

where the perturbation introduced by the non-minimal 
coupling is clearly visible. 

In order to scrutinize the profile of pressure and den- 
sity inside the Sun, one requires a suitable equation of 
state. Instead of pursuing a realistic representation of 
the various layers of the solar structure, one resorts to a 
very simplistic assumption, the so-called polytropic equa- 
tion of state. This is commonly given by p = JTp^""'"^^/"', 
where K is the polytropic constant, p is the mass density 
and n is the polytropic index. A polytropic equation of 
state with n = 3 was used by Eddington in his first solar 
model, and will be adopted here due. 

Given this equation of state, one may write p = Pcd^^{C) 
and p = Pc^"'''^(0) with ^ = t/tq a dimensionless vari- 
able and rl = {n+ l)pcl AnGpl; pc = 1.622 x 10'"^ kg/m^ 
is the central density, and Pc = 2.48 x 10^^ Pa is the cen- 
tral pressure. One obtains the perturbed Lane-Emden 
equation for the function 0{^): 



1 

e 



(41) 



C 3n- 1 



9' 4{n+l) 



4^ 2 
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FIG. 1: Relative perturbation 6 for 2.8 < n < 3.2. 



where the prime now denotes derivation with respect to 
the dimensionless radial coordinate ^, and one defines the 
dimensionless parameters A^. = Xpc and Nc = Pc/pc = 
1.7 X 10~^, for convenience. Clearly, setting = one 
recovers the unperturbed Lane-Emden equation [25] . 

Notice that the perturbed Lane-Emden equation is 
a third-degree differential equation; its numerical res- 
olution is computationally intensive and displays some 
complex behaviour; conveniently, the assumed perturba- 
tive regime prompts for the expansion of the function 
^(C) = ^^o(l + Ac6) around the unperturbed solution 
OoiO- Inserting this into Eq. (41) and expanding to 
first-order in Ac, one obtains 

S" + 2 f J + i V' + - 1)^0-''^ = ^-?^0r-% (42) 



+ (2n + l)Nc^0l--X + 4^i±402„-i 



2n 



5n(n- l) ^g„-3g/3 , "W" ' ; an-2a/2 



4(n+l) ° 

4(n+l) ^0 ^0 + 2~ 



1/1/ 
0*0 



supplemented by the initial conditions (5(0) = ^'(0) = 
0. Notice that the choice for the perturbative expansion 
leads to a solution 6 independent from the parameter Ac. 

After dealing with the issue of exterior matching con- 
ditions and bypassing a troublesome divergence of 5 near 
the boundary of the star [17], one may obtain the nu- 
merical solution for Eq. (42) for a polytropic index in 
the vicinity of n = 3, as depicted in Fig. 1. 

Finally, one turns to the issue of calculating one of 
the observables under scrutiny, that is, the central tem- 
perature of the Sun. The polytropic equation of state 
indicates that p oc T"+-^, which yields 



Tco 
Ac /■« 

Cr^Or Jo 



n+1 



(43) 



n6 ■ 



3n 9'^ 



2^ 8 



where = Rr/ro and = 0.713i?Q marks the onset of 
the convection zone (where the chosen equation of state 



FIG. 2: Relative deviation of the central temperature 
Tc/Tco — 1, with contour lines of step 0.1%. 



fails) and T^o is the central temperature derived from 
the Ac = unperturbed scenario. One may derive a 
parameter plot in the (n, Ac) parameter space, shown in 
Fig. 2. As can be seen, no relative deviation of the central 
temperature occurs above the experimentally determined 
level of 6%. However, since the values found are of the 
order of 1%. one may hope that any future refinement of 
the experimental error of Tc could yield a direct bound 
on the parameter Ac- Furthermore, the perturbative con- 
dition A Kpc is confirmed (reintroducing the factor k, 
for clarity), which translates to |A| < 4.24 x 10^^ eV"^. 



VI. MODELS WITH ARBITRARY COUPLINGS 
BETWEEN MATTER AND GEOMETRY 

The discussed gravity models with linear coupling be- 
tween matter and geometry, given by Eq. (1), can be 
further generalized by assuming that the supplementary 
coupling between matter and geometry takes place via 
an arbitrary function of the matter Lagrangian Lm, so 
that the action is given by [18] 

S = j^\h{R) + Q {Cm) [1 + A/2 (i?)]| s/^gd^x, 

(44) 

where Q {Cm) is an arbitrary function of the matter La- 
grangian density Cm- The action given by Eq. (44) 
represents the most general extension of the Einstein- 
Hilbert action for GR, S' = / [R/2 + Cm] ^f^d'^x. For 
fi{R) = i?, f2{R) = and Q {Cm) = Cm, one recovers 
GR. With f2{R) = and Q {Cm) = Cm one obtains the 
f{R) generalized gravity models. The case Q [Cm) = Cm 
corresponds to the linear coupling between matter and 
geometry, given by Eq. (1). The only requirement for /j, 
i = 1,2 and Q is that they are analytical functions of the 
Ricci scalar R and Cm, respectively - that is, they can be 
expressed as a Taylor series expansion about any point. 
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The field equations corresponding to action (44) are 

-2A {g^^n - V^V,) g (£„) F^iR) 

- [1 + \h{R)] [K - a {Cm)] g^. 

-[l + A/2(i?)]i^(/:„)V, (45) 

where F,{R) = df,{R)/dR, i = 1,2, and K {Crn) = 
dQ {Cm) /dCm, respectively. 

By taking the covariant divergence of Eq. 
(45), with the use of the mathematical identity 
[a'{R)R^, - a{R)g^,/2 + {g^,U - V^V.) a{R)] = 
[26], where a{R) is an arbitrary function of the Ricci 
scalar and a'{R) = da/dR, we obtain 

VT^. = Wln{[l + Xf2{R)]K{Cm)}{Cmg^,.-T^.) 

3C 

= 2Wnn{[l + Xh{R)]K{Cm)}^^. (46) 

For Q {Cm) = Cm, one recovers the equation of motion 
of massive test particles in the linear theory, Eq. (7). As 
a specific model of generalized gravity models with arbi- 
trary matter-geometry coupling, one considers the case in 
which the matter Lagrangian density is an arbitrary func- 
tion of the energy density of the matter p only, so that 
Cm — Cm (p) • One assumes that during the hydrody- 
namic evolution the energy density current is conserved, 
{pU") = 0. Then, the energy-momentum tensor of 
matter is given by 



T"^ = p^U^U-" +(Cm- P^ 
dp \ dp 



(47) 



where we have used the relation 5p = 
(1/2) p{gij^„ — U^Ui,)dg^'^, a direct consequence of 
the conservation of the energy density current. 

The energy-momentum tensor given by Eq. (47) can 
be written in a form similar to the perfect fluid case if 
one assumes that the thermodynamic pressure p obeys a 
barotropic equation of state, p = p{p)- In this case the 
matter Lagrangian density and the energy-momentum 
tensor can be written as 

Cm {p) = p[l + U {p)] =p(^l + J^'±^-p{p), (48) 
and 

T'^'' = {p [1 + n (p)] + p {p)} U^U^ +p{p) g^\ (49) 
respectively, where 



n(,)= r*-^. 

Jo P P 



(50) 



Physically, 11 (p) can be interpreted as the elastic (de- 
formation) potential energy of the body, and therefore 



Eq. (49) corresponds to the energy-momentum tensor of 
a compressible elastic isotropic system. From Eq. (46), 
one obtains the equation of motion of a test particle in 
the modified gravity model with the matter Lagrangian 
an arbitrary function of the energy density of matter as 
Eq. (6), where the extra force is now given by 



r = V.ln<^[l + A/2(i?)]i^[/:„ {p)\ 



dCm {p) 

dp 



(51) 

It is easy to see that the extra- force /'', generated due to 
the presence of the coupling between matter and geome- 
try, is perpendicular to the four-velocity, f^U/j, = 0. The 
equation of motion, Eq. (6), can be obtained from the 
variational principle 



5Sp = 5 J Cpds = 5 J ^/Q^/g^JMPds = 0, 



(52) 



where Sp and Cp 



Q\/ gtivU^^U" are the action and 
Lagrangian density, respectively, and 

^=[l + Xf2{R)]K[Cm{p)]^^^. (53) 

The variational principle Eq. (52) can be used to 
study the Newtonian limit of the model. In the 
weak gravitational field limit, ds « -y^l -t- 20 — iPdt « 
(l + — ir^/2) dt, where is the Newtonian potential 
and V is the usual tridimensional velocity of the parti- 
cle. By representing the function as 

v/Q = [l + Xh{R)]K[Cm{p)]^^^ 



= l + ^[R,Cm{p) 



dCm (p) 

dp 



(54) 



where |$| << 1, the equation of motion of a test particle 
can be obtained from the variational principle 



^[R,Cm (p), 



dCm (p) 

dp 



and is given by 



a = —V<j) — V$ = dN + CLE, 



dt = 0, (55) 



(56) 



where ojv = — V0 is the usual Newtonian gravitational 

acceleration and Se = — V$ a supplementary effect in- 
duced by the coupling between matter and geometry. 

An estimative of the effect of the extra-force generated 
by the coupling between matter and geometry on the or- 
bital parameters of planetary motion around the Sun can 
be obtained by using the properties of the Runge-Lenz 
vector, defined as ^ = vx L — acr, where v is the velocity 
relative to the Sun, with mass Mq, of a planet of mass m, 
r = rSr is the two-body position vector, p = pvis the rel- 
ative momentum, p = ttiMq/ (to + Mq) is the reduced 
mass, L = r X p = pr^Ok is the angular momentum. 
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and a = GitiMq [27]. For an elliptical orbit of eccen- 
tricity e, major semi-axis a, and period T, the equation 
of the orbit is given by (L^//xq;) = 1-1- ecos6'. The 
Runge-Lenz vector and its derivative can be expressed as 




fLee, 



and 



dA 



dV{r) 
dr 



ee, 



(57) 



(58) 



respectively, where V{r) is the potential of the central 
force [27]. The potential term consists of the Post- 
Newtonian potential, 



Vpjv(r) = 2' 



(59) 



plus the contribution from the general coupling between 
matter and geometry. Thus, one has 



dA 
Id 



+ maE{r) 



(60) 



where it is also assumed that fi fa m. The change in 
direction A<p of the perihelion for a variation of 9 of 2tt 
is obtained as 



1 



ae 











lo 





and is given by 



247r3 

2n as 



de, 



(61) 



I 



L"^ (1 + ecos^)~^ /ma 



{i + ecosey 



cos 9d9, 



(62) 



where the relation a/ L = 2t: {a/T) / \/\ — e? is used. The 
first term of this equation corresponds to the GR predic- 
tion for the precession of the perihelion of planets, while 
the second gives the contribution to the perihelion pre- 
cession due to the presence of the new coupling between 
matter and geometry. 

As an example of the application of Eq. (62), one 
considers the case for which the extra- force may be 
considered constant — an approximation that might be 
valid for small regions of the space-time. Thus, through 
Eq. (62), one obtains the perihelion precession 



67rGAf,5 2T^a^^fY^ 
a (1 - e^) GMq 

resorting to Kepler's third law, = Air'^a^ /GMq. 



(63) 



For Mercury, a = 57.91 x 10^ m and e = 0.205615, 
respectively, while Mq = 1.989 x 10'^'^ kg: the first term 
in Eq. (63) gives the GR value for the precession angle, 
(A^)g^ = 42.962 arcsec per century, while the observed 
value is (A0)^^^ = 43.11 ± 0.21 arcsec per century [28]. 
Therefore, the difference (Ai/))^ = (Ai/))^^^ — (A^)^,^ = 
0.17 arcsec per century can be attributed to other physi- 
cal effects. Hence, the observational constraints requires 
that the value of the constant extra acceleration ue must 
satisfy the condition 



aE < 1.28 X 10"^^ m./s. 



(64) 



This value of a_E, obtained from the solar system obser- 
vations, is somewhat smaller than the value of the extra- 
acceleration ao ~ 10~^° m/s^, necessary to account for 
the Pioneer anomaly [8]. However, it does not rule out 
the possibility of the presence of some extra gravitational 
effects acting at both solar system and galactic scale, 
since the assumption of a constant extra-force may not 
be correct on large astronomical scales. 



VII. CONCLUSIONS AND OUTLOOK 

In this contribution we have discussed a wide range of 
implications of the gravity model action, Eq. (1), whose 
main feature is the non- minimal coupling between curva- 
ture and the Lagrangian density of matter (or a function 
of it, in Section VI). This exhibits an extra force with re- 
spect to the GR motion, as well as the non-conservation 
of the matter energy-momentum tensor. The prevalence 
of these features for different choices for the matter La- 
grangian density was discussed in Section III. In Section 
IV, the specific features of the associated scalar-tensor 
theory were discussed — and it was shown that the model 
is consistent with the observational values of the PPN pa- 
rameters, namely /3 = 7 = 1, to zeroth-order in A. In Sec- 
tion V, we consider the impact of the novel coupling on 
the issue of stellar equilibrium. It is shown that, for the 
simplest model of the Sun, the effect of the new coupling 
on the central temperature is smaller than 1 %, which 
is consistent with the uncertainty of current estimates. 
Finally, in Section VI, a general function of the matter 
Lagrangian density has been introduced, and the value 
of the resulting extra force obtained, < 10^^^ m/s^. 

Of course, further work is still required in order to 
quantify the violation of the Equivalence Principle intro- 
duced by the model under realistic physical conditions. 
A low bound for the coupling A, would justify the results 
discussed in this work, which are first order in A. Impli- 
cations of the discussed model in what concerns the issue 
of singularities are still to be addressed, as well as the 
impact that the new coupling term might have on the 
early Universe cosmology. 

We would like to close this contribution with our best 
wishes to our colleague Sergei Odintsov, on the occasion 
of his 50th birthday. 
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